Observation of scalar nuclear spin-spin coupling in van der Waals molecules 
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Scalar couplings between covalently bound nuclear spins are a ubiquitous feature in nuclear mag- 
netic resonance (NMR) experiments, imparting valuable information to NMR spectra regarding 
molecular structure and conformation. Such couplings arise due to a second-order hyperfine in- 
teraction, and, in principle, the same mechanism should lead to scalar couplings between nuclear 
spins in unbound van der Waals complexes. Here, we report the first observation of scalar couplings 
between nuclei in van der Waals molecules. Our measurements are performed in a solution of hyper- 
polarized ^■^^Xe and pentane, using superconducting quantum interference devices to detect NMR 
in 10 mG fields, and are in good agreement with calculations based on density functional theory, 
van der Waals forces play an important role in many physical phenomena, and hence the techniques 
presented here may provide a new method for probing such interactions. 



Introduction 

Scalar coupling between nuclear spins, characterized 
by an interaction of the form Hj = hJli ■ I2, where h 
is Planck's constant, J is the coupling constant, and Ii 
and I2 are the two spins involved, often referred to as 
J-coupling, is a staple feature of NMR spectroscopy. fl| 
J-coupling arises primarily from a second-order Fermi- 
contact interaction, which is generally thought to require 
covalent bonding between the two coupled spins. In co- 
valent molecules, this interaction gives rise to observable 
splittings in NMR spectra whose pattern and separation 
provides important information on molecular structure 
and conformation. Such couplings vanish quickly upon 
increasing the number of intervening bonds, and cou- 
plings extending over more than 4-5 bonds are generally 
very weak if at all observable. As a consequence, the 
observation of a coupling constant is taken to be a sig- 
nature of covalent bonding, and is accordingly exploited 
to establish molecular connectivities. 

The two coupled nuclei are therefore generally as- 
sumed to belong to the same molecule. However, strictly 
there is no such requirement, as long as the electronic 
structure of the system allows for the two spins to be 
"connected" by a (possibly low) electron density. ^] Very 
long-range couplings involving ^^F have been known for 
a long time and often denoted as occurring "through 
space" (implying that the through-bond connectivity is 
too long to allow for coupling to be transmitted) . [1, 0] 
Hereafter, we will use the term "through space" even 
though the spin-spin interactions involved are the same 
as in covalent molecules. [2] Such a situation holds also for 
even relatively weak hydrogen bonds (HB), and through- 
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HB couplings are now routinely analyzed in structural 
biochemistry. [5] Similarly, through-space J{^^C,^^C) or 
J(^H,^H) couplings have recently been reported for pro- 
teins whose conformation maintains a sufficiently small 
distance between coupled nuclei [6] and for [2,2]paracy- 
clophanes Analysis of such couplings in molecules 
where the interacting spins are connected by covalent 
bonds require careful consideration of conformational ef- 
fects [1] and coupling pathways [oj in order to assess their 
nature. 

Scalar couplings in van der Waals molecules are ex- 
pected to arise through the same mechanism as in the 
case of covalently bound spins. Pioneering computational 
work on Xe- • • Xe and Xe- • • H predicted very small 
couplings (in the range of 10~^ Hz). Further works using 
improved theoretical methods, both ab initio and based 
on density functional theory (DFT) , have predicted cou- 
plings on the order of 10^^ Hz in He2, (CH4)2, a variety 
of organic [8, 11-14] and xenon [15] complexes. Similar 
calculations have enabled accurate prediction of through- 
HB couplings, fslliel 

Experimental detection of couplings in the range of 
10^^ Hz is feasible with conventional high-field NMR 
techniques. However, for unbound systems, chemical 
exchange and diffusion averages the effects of such cou- 
pling to very nearly zero in thermally polarized samples. 
Thus, despite the fundamental and practical implications 
of the existence of spin-spin couplings in unbound van der 
Waals molecules, there is no record of their experimental 
observation. 

In this work, we report the first observation of scalar J- 
couplings between unbound spins in van der Waals com- 
plexes in a solution of hyperpolarized liquid Xe and pen- 
tane. The approach is based on earlier observations [13] 
that J-coupling in the presence of fast chemical exchange 
can be observed if one spin species is hyperpolarized, re- 
sulting in an average frequency shift of the other species. 
With this technique, the strength of proton NMR signal 
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is enhanced by a factor of 10^ by SPINOE with hyperpo- 
larized ^^^Xe.[ll, [l^ The measured value of J-coupling 
is in good agreement with density functional calculations 
averaged over an ensemble of spins simulating the bulk 
liquid phase. Van der Waals forces play an important 
role in many physical phenomena, and hence the ability 
to detect scalar couplings due to such forces may prove 
extremely valuable. For example, since J-couplings in un- 
bound systems vanish quickly with distance, their detec- 
tion and understanding would provide a spatial restraint 
akin to the nuclear Overhauser effect. 



Results and Discussion 

Consider a system of two different spins, 1 (pentane 
protons) and 2 (^^^Xe nuclei). The Zeeman contribu- 
tion to the proton Hamiltonian in the presence of a mag- 
netic field B = is Hz — — ?i7ili • Bq, where 71 is 
the proton gyromagnetic ratio. In the presence of rapid 
chemical exchange, the normal J-coupling Hamiltonian, 
/iJIi-l2, averages to Hj = h{J)Ii{l2z), where (J) repre- 
sents the thermodynamically averaged coupling between 
^^^Xe and [TtI, and we have assumed that the xenon 
spin polarization is nearly parallel to the magnetic field 
(note that J is given in Hz here). The Larmor preces- 
sion frequency of protons (assumed to be positive) is thus 
'^i = (7i/27r)i?2 -|- Avi, where 

Al/l =-(J)(/2.) (1) 

and the minus sign originates from the difference in the 
signs of the Zeeman and J-coupling Hamiltonians. (J) is 
given by jTrj : 

(J) J Jir)exp{-V{r)/kT)(fr (2) 

where n2 is the number density of spin I2, J{r) is a 
distance-dependent coupling constant, and V{r) is an in- 
teratomic potential. Equation [5] can be rearranged so as 
to highlight its relationship with the bulk liquid struc- 
ture, characterized by a radial distribution function g{r): 

(J) =n2 I 4Trr^J{r)g{r)dr. (3) 
Jv 

The frequency shift Avi is analogous to that observed 
in optical pumping experiments; in vapor mixtures of 
alkali metals and noble gases, spin polarization of the 
unpaired electron of the metal produces a shift in the 
NMR resonance of the noble gas. The shift is due to the 
Fermi Contact interaction between the unpaired electron 
of the metal and the nucleus of the noble gas, mediated by 
the hyperfine coupling constant. [2^ For the case of two 
nuclear spins the frequency shift can be parameterized as 

Ai^i = ^n^M2^, (4) 
zvr 6 



where M2Z = {l2z) IJ'2n2 / 12 is the magnetization of 
spin 2 and /i2 is its magnetic moment (/12 = —3.88 x 
lO"^** erg/G in the case of ^^^Xe[2l|). k is a dimen- 
sionless parameter representing the suppression or en- 
hancement of the classical dipole magnetic field due to 
Coulomb repulsion or attraction of the particles. Recall- 
ing Eq. dl]) and setting I2 ~ 1/2, we have 

g 

(J) = --K7i/i2n2. (5) 

SQUID Measurements. Our measurements are per- 
formed using hyperpolarized xenon (typical polarization 
is 2-3%). Despite the use of hyperpolarized ^^^Xe, fre- 
quency shifts of Eq. [T] are quite small, on the order 
of tens of mHz, making detection in a high field spec- 
trometer challenging, due to the finite linewidth of res- 
onances (either because of homogeneous or inhomoge- 
neous broadening) and due to small magnetic field drifts. 
We work in a low-field (10 mC), magnetically shielded 
environment (see Supplementary Information, Fig. 1), 
where the high absolute field stability facilitates obser- 
vation of such small frequency shifts. NMR signals are 
monitored via high Tc superconducting quantum inter- 
ference devices (SQUIDs) to achieve high sensitivity at 
low frequencies (magnctometric sensitivity was roughly 
0.3 nG/\/Hz). Our measurements were performed us- 
ing the following protocol: 1) The spherical sample cell 
was filled half full of pentane, which was then degassed 
via three freeze-thaw cycles under vacuum, and finally 
held in liquid state at about — 100°C. 2) Hyperpolar- 
ized xenon gas was introduced into the sample cell where 
it condensed to liquid state to fill the remaining half of 
the cell. 3) Spin precession was excited by applying a 
7r/2 pulse to the protons via a rotating magnetic field 
at the proton resonance frequency. Although this pulse 
was off resonance for the ^^^Xe spins, it produced small 
excitation , tipping them into the transverse plane by 
about 1°. Proton T2 was typically about 3.2 s, while 
^^^Xe Ti was about 450 s, so that many (« 50) proton 
transients could be collected for each batch of polarized 
xenon. Free induction decay signals for ^H and ^^^Xe nu- 
clei are shown in Fig. [Ija), after passing the raw signal 
through bandpass filters centered about the Larmor pre- 
cession frequencies at 39.2 and 10.8 Hz, respectively. The 
magnitude Fourier transform of the raw signal is shown 
in Fig. [TJb) where the ^^^Xe and proton resonances can 
be easily identified. 

The signature of the desired interaction is a proton 
frequency shift that is linear in the xenon magnetiza- 
tion. There are two sources of spurious frequency shifts 
that need to be accounted for. We first address dipo- 
lar fields associated with a deviation of the cell from 
perfectly spherical geometry (primarily due to the filling 
port). These effects can be compensated for by acquiring 
data with the magnetic field oriented in three orthogonal 
directions. We acquired three data sets with Bq oriented 
along each of the three directions y, and z. 

Figure m a) shows the decay of the longitudinal compo- 
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FIG. 1: Typical signals obtained in our experiment after ap- 
plying a 7r/2 pulse to the protons, (a) Free induction decay 
signals for ^^^Xe and after passing the raw signal through 
the appropriate bandpass filters. Proton T2 was 3.2 s for these 
data, (b) Magnitude Fourier transform of the raw signal. Ow- 
ing to the very small applied field, proton chemical shifts in 
pentane are unresolved, yielding a single proton line at 39.2 
Hz. 



nent of the magnetization, determined from the DC value 
of the SQUID signals. Overlying these data is a decaying 
exponential with time constant Ti k, 450 s. -'^^^Xe and 
proton resonance frequencies following each proton 7r/2 
pulse were determined by passing the raw signal through 
the appropriate bandpass filter and fitting the result to 
a decaying sinusoid. Figure [2][b) shows the deviation of 
the and ^^^Xe NMR frequencies from their value in 
the limit of low ^■^^Xe magnetization as a function of the 
^^^Xe magnetization. Shown here are typical data sets 
for Bo oriented in three different directions. 

The second source of spurious frequency shifts is due to 
the magnetic fields associated with the pickup and feed- 
back coils of the SQUID magnetometers. The size of this 
effect can be accurately determined based on the geome- 
try of the detectors and several auxiliary spin-precession 
measurements (see Methods and Supplementary Infor- 
mation). After correcting for the effect of the SQUIDs, 
we verify that -'^^^Xe frequency shifts vanish (scalar spin- 
spin couplings between ^^^Xe spins are not observable), 
and extract the part of the proton frequency shift that is 
attributed to scalar coupling with ^^^Xe. 

Frequency shifts are summarized in Table U the last 
row gives the net frequency shift after subtracting the 
shifts due to the SQUID magnetometers. The net ^^^Xe 
frequency shift is consistent with zero, and the average 
proton frequency shift per unit of magnetization along 
the three axes is Aj/i/M2z = — 52±11 Hz/G, correspond- 



FIG. 2; (a) Decay of longitudinal component of xenon mag- 
netization, monitored by the DC component of the SQUID 
signal in a typical data set. (b) Frequency shifts of ^^®Xe and 
NMR signals for the magnetic field oriented in three or- 
thogonal directions. The straight lines overlying the data are 
linear fits. 



^^^Xe, corresponding to a density n2 
Employing Eq. ([5]), we find (J) = 



ing to K = -0.0014± 0.0003. The cell was half-filled with 

7x10^1 cm-3.[22| 
2.7 ± 0.6 Hz. We 
note that accuracy and precision can be considerably im- 
proved by (a) an apparatus that overcomes some techni- 
cal limitations (cell distortions, variations in cell geom- 
etry due to filling, and inaccurate measurement of the 
longitudinal magnetization for certain configurations of 
the magnetic field); (b) low- Tc SQUIDs with better sen- 
sitivity (« 0.01 nG/v^) than the high-T^ SQUIDs used 
here; (c) a higher xenon polarization than the 2-3% used 
in this work (polarization as high as 65% has been re- 
ported [H). 

The measured value of (J) cannot be directly com- 
pared with DFT-calculated coupling constants for a sin- 
gle Xe-pentane van der Waals complex, since in the bulk 
phase a variety of molecular arrangements exists; thus 
we have combined DFT calculations with a simulation of 
the liquid phase as follows. 

Computational studies. Two separate sets of calcu- 
lations were carried out: First, through-space couplings 
were calculated with relativistic DFT methods for a se- 
ries of xenon-pentane arrangements, thereby character- 
izing the most relevant portion of the ^^^Xe-^H coupling 
surface, see Figure 2 of Supplementary Information. Sec- 
ond, we determined the radial distribution function g{r) 
for the bulk liquid phase, which allows averaging of cal- 
culated coupling constants via Eq. Most calculated 
couplings are negative and lie between —8 and Hz (a 
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Bo direction 




Aiy2/M2z 


Aiy2,pure/M2z 


X 


-160 ± 10 


-21 ± 1 


-13 


y 


162 ± 29 


106 ± 12 


96 


z 


-170 ±9 


-32 ±2 


-31 


x+y+z 


-168 ±32 


53± 12 


53 


SQUID shift 


-12 


39 


39 


net 


-156 ±32 


14± 12 


13 



TABLE I: List of frequency shifts per unit of magnetization 
(Hz/G) for the magnetic field oriented in orthogonal direc- 
tions. Errors are determined from the scatter of points over 
several trials, and are not given for pure xenon runs because 
we have collected only one data set for each direction. 



few positive values are also found, as previously reported 

The magnitude of these couplings is consistent with the 
results of other computational works based on various ab 
initio and DFT theoretical methods jsl. [lll-[T5| . However, 
our values are larger than those reported by Salsbury and 
Harris for Xe- • • H 10] by 2-3 orders of magnitude. In this 
regard, we note that this seminal paper employed 
the Thonias-Fermi-Dirac functional, which is based on 
the uniform electron gas model. This early functional is 
known to suffer from neglect of electron correlation and 
inaccurate description of the electron density found in 
molecular systems [24j . 

The main trend of the calculated couplings features a 
roughly exponential decay as the distance increases from 
ca. 3 to 4 A (See Supplementary Information, Figure 3). 
Accordingly, the couplings from the second solvation shell 
are mostly negligible. An orientational dependence of the 
coupling constants is evident from the spread of points: 
some relatively large couplings are predicted even when 
the Xe-H separation is larger than 4.0 A. Similar anoma- 
lous distance dependences were observed in Xe- • • CH4 
[is"] and in HF- • • CH4 [ll|, ^ complexes. A detailed ex- 
ploration of the angular dependence would probably not 
provide important new information, since we are only 
interested in the average properties of the bulk liquid 
phase. Calculated couplings in the range 3.0-4.3 A were 
fitted to an exponential decay curve J(r) — ae^'"', with 
a = -2.7867 x lO"* Hz, b = 2.719 A-\ 

In Figure [3] we report the radial distribution function, 
g(r), of the distance between xenon and the hydrogen 
atoms of pentane together with the calculated couplings 
and the fitting curve — J(r). There are three groups 
of chemically inequivalent protons in pentane, however 
since chemical shifts are unresolved in our experiment, 
we average the radial distribution function over all hy- 
drogen atoms. As we can see, the probability to find a 
Xe-H pair of atoms separated by less than 3 A is mostly 
negligible, thus there is no need to carefully characterize 
the J{r) function for such short distances. Integration 
of the function of Eq. Hwith = 216/(31.76 A)^ (see 
Simulation Details below and in Supplementary Infor- 
mation) gives an average value (J) = —3.2 Hz, within 




3.5 

R / Angstrom 

FIG. 3: (Thin line) Radial distribution function, g(r), of the 
Xe-H distance averaged over all protons of pentane, right axis. 
(Open circles) calculated coupling constants with a negative 
sign for easier comparison with g{r). (Thick line) fitting curve 
-J(r) = ae-'"', with a = -2.7867 x 10* Hz, b = 2.719 A-\ 
left axis. 



experimental error of SQUID measurements. 

It is worthwhile to stress that the averaged coupling 
that we have determined is somewhat different from that 
encountered in high-field NMR spectra, where it is usu- 
ally possible to measure the coupling constants between 
individual spin pairs. In contrast, the coupling we have 
determined between pentane protons and xenon is the 
sum of all coupling constants between protons and all 
xenon atoms surrounding pentane, since each coupling 
produces a shift which adds to the shifts of the others. 
This situation is illustrated by the MD snapshot in Figure 
m where we have highlighted a single pentane molecule 
and the xenon atoms in the first solvation shell. Given 
the rapid decrease of J(r) with distance, only the closest 
xenon shell provides a sizable contribution to (J). Thus, 
the value of (J) is close to that calculated at a Xe- • -H 
distance of ca. 3.2 A. 

For the same reasons, (J) depends on the density. Al- 
though the radial distribution function profile may also 
be density-dependent, the largest effect is simply a linear 
scaling of the coupling with the density of hyperpolarized 
xenon. Thus, in the limit of infinite dilution of pentane 
(cell filled with xenon), the frequency shift, and there- 
fore (J), would be twice as large as in our case, where 
the mixture is only 50% xenon. 



Conclusions 

We have reported the observation of a scalar coupling 
between ^^^Xe and ^H in a solution of hyperpolarized 
xenon and pentane. The measured value is in good 
agreement with calculations based on density functional 
theory, averaged according to the bulk liquid structure 
obtained by simulation. To the best of our knowledge. 
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Methods 




FIG. 4: Snapshot of the MD simulation box where we have 
highlighted a single pentane molecule and the closest xenon 
atoms more strongly contributing to the measured frequency 
shift of the pentane protons. 



this is the first observation of scalar coupling mediated 
by van der Waals interactions in the presence of rapid 
chemical exchange, i.e. for "unbound" spins. Given the 
rapid decay of such couplings with the xenon-probe dis- 
tance, their observation provides cogent information on 
the spatial proximity of unbound partners without com- 
plications arising from relaxation effects, as is often the 
case with nuclear Overhauser effects. Our measurements 
are facilitated by operating in very low magnetic fields, 
allowing us to easily vary the direction of the magnetic 
field Bq over three orthogonal orientations to average 
away the effects of dipolar fields. In principle magic angle 
spinning in high field may accomplish the same thing, al- 
though the necessity of vacuum tight mechanical connec- 
tions to the xenon hyperpolarization system would make 
this challenging. Implementation of pulse sequences such 
as WAHUHA (Waugh-Huber-Haberlein) ^ , designed to 
mitigate the effects of dipolar fields may also enable high- 
field detection of such effects. 

M.P.L., N.T. and M.V.R. designed and carried out 
SQUID experiments, analyzed the data and wrote the 
paper. A.B. and G.S. designed and carried out DFT 
calculations and MD simulations, analyzed the data and 
wrote the paper. 



Experimental SQUID measurements. As typi- 
cal for DC SQUID magnetometers [2d\, ours consist of 
a SQUID loop with two Josephson junctions, a super- 
conducting pick-up coil and a normal feedback coil (inset 
(a). Figure SI in the Supplementary Information). A 
feedback circuit controls the current in the feedback coil 
such that the SQUID is locked to a fixed point of its 
V{^) curve. The current in the feedback coil serves as 
a measure of the applied magnetic field, but also gener- 
ates small magnetic fields which can affect the precession 
frequency of the nuclear spins. We consider separately 
the AC and DC magnetic fields generated by the SQUID 
coils. The AC field is generated in response to the oscil- 
lating z component of the magnetization (inset (b). Fig. 
SI in the Supplementary Information), and is in-phase 
with Mz- It can be shown from Bloch equations that the 
AC in-phase field will cause a shift of the ^^^Xe NMR 
frequency by slowly rotating the longitudinal component 
(parallel to Bq) of the Xe magnetization into the trans- 
verse plane. However, it does not cause an appreciable 
shift of the proton NMR frequency because the proton 
spins are tipped by approximately tt/2 relative to Bq. 
In contrast, the DC field will affect both spin precession 
frequencies. It is important to take into account the fact 
that the SQUID feedback circuit is reset after every NMR 
pulse, which results in a much smaller DC field. We di- 
rectly measured the AC and DC fields of the SQUIDs 
at the location of the spin sample by carefully measur- 
ing the dependence of the ^^^Xe NMR frequency on the 
applied uniform magnetic field. 



DFT calculations and molecular dynamics sim- 
ulation. Spin-spin coupling involving heavy-atom nu- 
clei such^ as xenon are affected by major relativistic 
effects, [27| which have been dealt with by the imple- 
mentation of the Zero-Order Regular Approximation 
(ZORA) method in the framework of DFT using the soft- 
ware package ADF.[28| Thus, we ran a series of calcu- 
lations of the through-space J(^^^Xe,^H) coupling con- 
stants between the protons in a pentane molecule and a 
xenon atom, at the scalar ZORA BP86/TZ2P level, con- 
sidering only the Fermi contact, diamagnetic spin-orbit, 
and paramagnetic spin-orbit terms, pol l30l| The position 
of the xenon atom was varied over a cubic grid of points, 
with the origin at the center of mass of the pentane 
molecule for a total of 1428 calculated coupling constants 
(see Supplementary Information, Fig. 3). The radial dis- 
tribution function in the bulk liquid phase was obtained 
by molecular dynamics simulations, for which we used 
an all-atom force field (OPLS-AA) for pentane [3l| and 
a Lennard- Jones model for xenon [s^ (see the Supple- 
mentary Information for details). 
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